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Abstract 

We recall our recent description of quark parton densities of the proton 
at Q 2 = 4 GeV 2 in terms of Fermi-Dirac distributions parametrized with 
very few free parameters. We have also proposed some simple assumptions 
to relate unpolarized and polarized quark parton densities which lead to 
a fair description of the spin-dependent structure functions xgi(x,Q 2 ) and 
xg™(x, Q 2 ) at low Q 2 . We will show the predictions we obtain after a straight¬ 
forward DGLAP Q 2 evolution and comparison in a much broader x and Q 2 
range, with several recent and accurate deep-inelastic scattering data. In 
particular, we will see that we get an excellent agreement with the sharp 
rise of F 2 ep (x, Q 2 ) for small x, recently observed at HERA. Finally, we give 
several predictions for lepton pair and gauge boson production in pp and pn 
collisions at high energies which will be tested in the future at RHIC. 
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1 Introduction 


Deep-inclastic lepton-hadron scattering is the basic source of our knowledge 
on quark parton densities of the proton and an enormous amount of data 
has been accumulated over the last twenty years or so. It allows to test 
rather accurately our picture of the nucleon structure at short distances. 
With the advent of HERA, a new kinematical range is now accessible in x 
down to 10 -4 and in Q 2 up to 10 4 GeV 2 , which will be extremely useful for 
the physics analysis at future hadron colliders. Besides the determination, 
for each flavor u,d,s, etc..., of the valence quarks which dominate, say for 
x > 0.1, and of the sea quarks (or antiquarks) which dominate in the small 
x region, it is very important to extract the correct gluon distribution. This 
is because, first it also prevails at low x and second, it plays a crucial role in 
the Q 2 evolution for testing perturbative QCD. 

Let us now briefly review the main points of our approach. Since quarks 
carry a spin—1/2, it is natural to consider that the basic distributions are 
xqf(x,Q 2 ) corresponding to a quark of flavor i and helicity parallel or an- 
tiparallel to the proton helicity. Recently we have proposed || a simple 
description of these quark densities in terms of Fermi-Dirac distributions 
with very few free parameters, which were determined from the data at 
Q 2 = Qq = 3 GeV 2 . This statistical physical picture for the nucleon struc¬ 
ture functions is largely motivated by the importance of the Pauli exclusion 
principle, which can be advocated to explain several experimental features, 
as discussed in ref. [jlj. The fact that there is also some experimental evidence 
for the existence of simple relations between unpolarized and polarized quark 
parton densities has allowed us to reduce to four, the total number of inde¬ 
pendent distributions, i.e. two for valence quarks xu± al (x,Ql) and two for 
sea quarks (or antiquarks) xuf ea (x,Q jj). In this approach, the corresponding 
four d quark densities xdff al {x, Qq) and xdf ea (x,Ql) are related to the four 
u quark densities mentioned above and similarly for the s quark (or anti¬ 
quark) densities. Concerning the gluon distribution xG(x, Q 2 ), for the sake 
of consistency, at Q 2 = Q q we have used a Bose-Einstein expression. As we 
will see, on the one hand, it is fully consistent with our present knowledge 
of xG(x, Q 2 ) extracted from deep-inclastic scattering data and, on the other 
hand, it leads to the correct Q 2 evolution of the structure functions. Actu¬ 
ally, one of the main purpose of this paper is to show that, starting from our 
simple scheme, which is in very good agreement with all spin-average and 
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spin-dependent structure functions for Q 2 near Q 2 = Ql, one gets, after a 
standard QCD Q 2 evolution, predictions fairly consistent with the data, in a 
much broader x and Q 2 range. 

The paper is organized as follows. In section 2, we briefly review and 
update our parametrization of the quark (antiquark) parton and gluon den¬ 
sities at Q 2 = Q o and show their Q 2 evolution. In section 3, the predictions 
we obtain for the spin-average structure functions F- 2 P (x,Q 2 ) with i = e,p 
and xF^ N (x,Q 2 ) are compared with various deep-inelastic scattering data, 
including the recent results from HERA. We also give our predictions for the 
W production cross section at the Tevatron and for the asymmetries in Drell- 
Yan lepton pair production and W production in pp and pn collisions, which 
are sensitive to the flavor asymmetry of light sea quarks in the proton. Sec¬ 
tion 4 is devoted to the discussion of quark and gluon helicity distributions 
at Q 2 = Qq and the comparison of our predictions with the most recent data 
on the spin-dependent structure functions xg p (x,Q 2 ) for proton, xg™{x, Q 2 ) 
for neutron and xgf(x, Q 2 ) for deuteron. We also discuss their Q 2 evolution 
and give our predictions for single and double helicity asymmetries for pp 
collisions at RHIC. Finally, in section 5 we consider the transversity distribu¬ 
tion of quarks (antiquarks) denoted by h\(x), which can be extracted from 
Drcll-Yan lepton pair and Z production in pp collisions with both proton 
beams transversely polarized. We give our concluding remarks in section 6. 


2 Quark (antiquark) parton and gluon distri¬ 
butions 


Let us first consider the distributions at Q 2 = Qq — 3 GeV 2 . For valence 
quarks we have two basic densities xu± al (x, Q jj) which are expressed in terms 
of Fermi-Dirac distributions of the form 


xp(x : Ql) 


n 

exp[(x — Xp)/x) + 1 


( 1 ) 


where Xp plays the role of the ’’thermodynamical potential”, which is a con¬ 
stant for each quark species and x is the ’’temperature” which is the same 
for quarks, antiquarks and gluons. This universal constant is independent of 
the parton helicity. We use a fitting procedure slightly different from that of 
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ref.[]I[ which is based on the most accurate neutrino data from CCFR [|. Q 
for xF^ N (x, Q 2 ), which has pure valence contributions, and also gives the 
antiquark distribution at Q 2 = 3GeV 2 (see below). Of course we will check 
the good agreement with the most recent NMC data|f| for F^(x,Q 2 ) and 
F?( x , Q 2 ) at Q 2 = AGeV 2 . It yields the following universal temperature and 
four free parameters entering in eq.(|l]) 

x = 0.092, b + = 0.354, b_ = 0.738, x*+ = 0.510, xl = 0.231 . (2) 

Clearly one has 

XUval ( x,Ql)= xu+ al (x,Ql)+ xu~ al {x, Ql) , (3) 

and, as we argued in ref. []TJ , we take 

xd va i(x, Ql) = 2 xu~ al (x, Ql) . (4) 

a± are not free parameters, but two normalization constants for the valence 
quarks in the proton. We note that the values of these parameters (see 
eq-fll)) are different from those of ref. (TH and correspond to a better fit of 
xF£ N (x, Ql) at large x, as shown in Fig.la. In particular we now find that 
the temperature x is slightly smaller than what we found in ref. [l]]. We also 
note that the potentials 5Q and xl for the two different hclicity states have 
numerical values in the ratio y/5 which is smaller than what we found in our 
previous analysis. As expected xQ > xl, so u+ al dominates over u! al . 

We now turn to antiquarks (or sea quarks) and for xh ± (x, Ql) we use the 
same expression eq. but in this case the potential has a smooth x depen¬ 
dence. This might reflect the fact that in this statistical description of the 
proton, we must consider the existence of two phases : a gaz corresponding to 
valence quarks, which dominates at large x, with a constant potential and a 
liquid corresponding to sea quarks (or antiquarks), which prevails at small x, 
with a potential slowly varying in x, that we take linear in y/x. In addition, 
we expect quarks and antiquarks to have opposite potentials, consequently 
the gluon which produces qq pairs has a zero potential (see below). More¬ 
over since in the process G —> q sea + q, q sea and q have opposite helicities, 
the potentials for uf ea (or u + ) and u~ (or uj ea ) will be opposite. So we take 

X U + = -Xu- = X 0 + Xl y/x . (5) 
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For the d distributions, we assume no polarization at Q q so we take 


xd + (x, Ql) = xd (. x,Ql)—xu (x,Ql) , (6) 


in accordance to eq.([|). Since when x —> 0, from Pomeron universality, one 
expects xu(x, Ql) = xd(x, Ql) ^ 0, a_ is not a free parameter and will be 
fixed by this constraint. This implies in addition b + = b_. For the strange 
quark (or antiquark) distributions, we also assume they are unpolarized and 
we take, in agreement with neutrino deep-inelastic scattering data, 

xs(x, Ql) = xs(x, Ql) = J (xu(x, Ql) + xd(x, Ql)) . (7) 

Therefore the antiquarks depend on four free parameters which are different 
from those determined in ref.Jl|, that is 


a + = 0.0185, b + = 6_ = —0.340, x$ = 0.219 and x\ = —0.406 . (8) 

We show in Fig. lb the result of our fit for xq(x) at Q 2 = 3GeV 2 . The 
NMC data for Ff and Ff which involve both valence quarks and sea quarks 
have also been used to test the correct determination of our four basic quark 
(antiquark) densities and the results are shown in Figs.2a,b. 

For comparison, we have also given in Figs.l and 2 the curves obtained 
by using the MRS(A) parametrization at Q 2 = AGeV 2 , taken from ref.||. 

Finally concerning the gluon distribution, for the sake of consistency, we 
use a Bose-Einstein expression given by 


xG(x,Ql) 


acx bG 
g x/x _ 


( 9 ) 


with a vanishing potential and the same temperature x as we discussed above. 
It is also reasonable to assume that for very small x, xG(x, Ql) has the same 
dependence as xq(x, Ql), so we will take bo = 1 + b, where b = b± is given 
in eq.(|8|). So except for the overall normalization o,q, xG(x,Ql) has no free 
parameter. From the momentum sum rule, we find ac = 13.146. For the 
sake of completeness, we also need to specify the gluon polarization for which 
there is no data at all, but we will make some simple speculations in section 4. 

To summarize, this statistical approach of the nucleon allows the con¬ 
struction of all quark, antiquark and gluon distributions in terms of simple 
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expressions which depend on nine free parameters. We gave some compre¬ 
hensive arguments about the physical meaning of these parameters for which 
we don’t have yet a full understanding. This small number of free parame¬ 
ters has been determined from deep-inelastic scattering data [Q, |3j. f|] at low 
Q 2 . Then, by means of a straightforward DGLAP || Q 2 evolution, we will 
compare our predictions with the existing data in a broad kinematical range 
of x and Q 2 , in order to test scaling violations and the dynamics implied by 
perturbative QCD. This approach contrasts with other methods presented 
in the literature f>j. 0 where one performs a global analysis with parton 
densities expressed in terms of a large number of free parameters, of the 
order of twenty, which are fixed by fitting several hundreds of data points. 
Notice that they are only dealing with spin-average distributions and it is 
necessary to introduce a new set of parameters to describe spin-dependent 
structure functions ||. Before closing this section, we would like to show the 
x-shapes of various unpolarized parton distributions at Q 2 = Ql = 3 GeV 2 
and Q 2 = 20 GeV 2 , after a standard Q 2 evolution. As already mentioned in 
ref.[]I[, we have used a numerical solution || of the DGLAP equations. We 
display in Fig.3 the u and d (valence + sea) quark distributions and also 
the corresponding antiquark distributions at two different values of Q 2 . The 
gluon density xG(x, Q 2 ) is shown in Fig.4. In the next section we will com¬ 
pare the results of our calculations, using these parton densities, with several 
pieces of existing experimental data and we will also give our predictions for 
future measurements. 


3 Experimental tests for unpolarized parton 
densities 

Since all parton densities have been determined, in order to test our approach, 
we can now proceed and calculate several physical quantities measured either 
in deep-inelastic scattering or in hadronic collisions. We will first study vari¬ 
ous structure functions in deep-inelastic scattering and then turn to hadronic 
processes with a special emphasize for testing the flavour asymmetry of the 
light sea quarks, i.e. u ^ d. 
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3.1 Deep-inelastic scattering 


We first consider the high statistics uN deep-inelastic scattering CCFR data 
p] from which one extracts xFg N (x, Q 2 ). As noticed above, this structure 
function gives a precise measurement of the isolated valence quark contribu¬ 
tions and we show in Fig.5 the results of our calculations. These neutrino 
data were obtained on a iron target, but we don’t think we can treat properly 
the nuclear effects. Therefore we have not tried to make any heavy nuclear 
target corrections which are known to be less important for xF£ N than for 
xFf N which contains the sea quark distributions. However in order to get 
the remarkable agreement displayed in Fig.5 we had to shift up by 4% our 
theoretical predictions. 

Next we turn to pp and ep deep-inelastic scattering for which several 
experiments have yielded a large number of data points on the structure 
function Fx,Q 2 ). First we will analyze the fixed target measurements 
which cover the limited kinematical region 0.0125 < x < 0.75 and 3 < Q 2 < 
230 GeV 2 , one obtains by combining the NMC data [10] and the BCDMS data 
ITT] . Our predictions are now compared with the data in Figs.6a,b but in 
this case we need to shift up by 12% the theoretical curves. The description 
of the data is rather impressive except at x = 0.75 for Q 2 < 10 2 GeV 2 where 
the BCDMS data he above the prediction. This is surprizing to us, because 
this kinematic region is dominated by valence quark contributions which are 
in perfect agreement with the measurement of xFf N at large x, as we have 
seen in Fig.5. 

We now turn to the very recent measurements of Ff p at Hera from the 
Zeus Collaboration 0 and the HI Collaboration Jl3| . It is important to note 
that these data are essentially in the low x region which is dominated by the 
sea quark densities. The behaviour of Ff p is therefore mainly constrained 
by the parameters given in eq.(^) and specially by b = —0.340, which was 
fixed by the steep behaviour of xq(x) at low Q 2 , as shown in Fig.lb. So 
this rise of the Hera data in the small x region was predictable from the 
CCFR measurement of xq(x) and no new parameter is requiered. This is 
not what was claimed in ref.|| where A = —0.3 was introduced at posteriori. 
Our predictions, with no overall normalization factor, are compared with the 
data in Figs.7a, b and c and the agreement is absolutely remarkable. From 
these new measurements of Ff p , by analyzing the scaling violations, it is 
possible to improve our knowledge on the gluon density in the low x region. 
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Going back to Fig.4, we see that xG(x, Q 2 ) increases very rapidly, and as we 
expect, it is fairly consistent with the distribution extracted at Q 2 = 20 GeV 2 
from the Zeus Collaboration [1141 . 


So far we have tested our parton densities in a kinematical region with 
Q 2 > Qq = 3 GeV 2 and we have seen that one gets a very satisfactory de¬ 
scription of all deep-inelastic structure functions. Clearly this confirms the 
success of the DGLAP Q 2 evolution for values above our starting point at 
Q 2 = Qq. Without trying to make a systematic study of the low Q 2 , low 
x region in electroproduction, we would like to show what we obtain for 
Q 2 < Qq, corresponding also to rather small x values. In this region there 


are already some data from NMC [il|, also earlier from SLAC [|15|] and some 
preliminary results from the FNAL muon scattering E665 experiment 
Just for illustration, we compare in Figs.8a and b our predictions with the 
NMC and SLAC data for Q 2 > 1 GeV 2 and our curves for x < 0.007 are 
given to be confronted with the final E665 data. For Q 2 < 1 GeV 2 the under¬ 
standing of the electroproduction structure functions certainly lies outside a 
perturbative QCD framework and involves other dynamical concepts JT7 . 


3.2 Hadronic processes 

As we discussed in ref. [|1] , one of the main features of our approach is the 
flavor symmetry breaking of the light sea quarks which was first recognized 
in ref.[[H| following an earlier NMC measurement [|TUJ . Here we will stress 
again the importance of this fact, we will indicate how one can confirm it 
and we will see where to expect ii(x,Q 2 ) ^ d(x,Q 2 ), for different x and Q 2 
values. 

Let us first consider the W production cross section which provides an im¬ 
portant test of the quark densities determined from deep-inelastic scattering. 
For illustration we show in Fig.9 the rapidity distribution for pp W + + X 
at Tevatron energy, which has been computed in the Drcll-Yan picture, dom¬ 
inated by the product u(x, M^)d(x, M 2 V ) with x ~ 0.05. We find that the 
maximum value of da/dyw is around 2.6 nb which is somehow smaller com¬ 
pared to other predictions presented in ref. . For the total W cross section 
we find % = 19 nb which is compatible with recent experimental values from 
CDF and D 0 at FNAL |29|. Next we turn to hadronic processes dominated 
by quark-antiquark annihilation, in order to study further the flavor symme¬ 
try breaking in the light sea quarks of the proton. One possibility, as noticed 
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in ref. [^lj , is to compare dilepton production in pp and pn collisions by means 
of the Drcll-Yan asymmetry 


A D y = 


d&pp/dy d(jp n ldy 
da pp I dy + da pn /dy 


( 10 ) 


These cross sections can be easily written in terms of q(x a , Ad' 2 ) and q(x b , Ad 2 ) 
and at rapidity y — 0, one has x a = xi, = \Jr = M/a/s, where M is the lepton 
pair mass. For the sake of simplicity, by neglecting the sea-sea contributions, 
one gets 

a _ (4A„ ~ 1)(A S ~ 1) + (A^ - 1)(4A S - 1) , . 

Dy (4A„ + 1)(A S + 1) + (A„ + 1)(4A S + 1) ’ 1 J 

where A„ = u va i/d va i and A s = u/d. In a broad kinematical region where 
A„ ~ 2, this expression shows the sensitivity of A DY to A s and, for example, 
for a flavor symmetric sea i.e. A s = 1, one gets A DY = 0.09, whereas A s = 0.5 
leads to Any = —0.09. This asymmetry has been measured at CERN by 


the NA51 Collaboration |22| in proton-proton and proton-deuteron collisions 
with 450GeI//c incident protons, in the mass range 4.3 < Ad < 8.5 GeV. As 
shown in Fig.10, the asymmetry Any predicted by our parton densities is 
compatible with the result of this rather low statistics experiment. A more 
accurate determination of A s should be obtained in the near future at FNAL 
by the E866 experiment [[£|. On the other hand, since there is a realistic 


possibility of having proton-proton and proton-deuteron collisions at RHIC 
with a high luminosity [^4], |23j, the measurement of Any should be seriously 
envisaged. In view of these future data, we give in Fig.11 the ratio (d — u)/ 
(d + u) versus x for a standard lepton pair mass Ad = 5 GeV (solid line). 

A second possibility for testing the value of A s has been proposed in 
ref. p5| by means of the production of W ± and Z in pp and pn collisions, also 
accessible at RHIC which will have the requiered energy. For W ± production, 
let us consider the following ratio 


R 


w 


_ dgW + /dy + dgW /dy - da^ + /dy - da™ /dy 
daZ + /dy + daZ~/dy + daZ + /dy + daZ~ /dy ' 


( 12 ) 


If we neglect the sea-sea contributions, it simply reads in terms of A v (x) and 
A s (a?) 

R _ (K{Xg) - l)(X s (x b ) - 1 ) + (x a ^ X b ) 

w (A„(x a ) + l)(A s (a b ) + 1) + (x a ^ x b ) 
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(13) 















Clearly Rw is symmetric under y —> —y and Rw = 0 if the sea is flavor 
symmetric, i.e. X s (x) = 1. We show in Fig.12 onr predictions for Rw at two 
different energies. We find that Rw is positive because A. s is always less than 
one, as we can check from the dashed line in Fig.11 which corresponds to 
Q 2 = M' 2 V . Rw decreases when y/s increases because for smaller x, according 
to Fig. 11, X s (x) increases and it follows from eq.([L3]) that Rw gets smaller. 
Actually the trend shown in Fig. 11 is simply due to the fact that the difference 
d—u does not change with Q 2 , whereas the sum d+u increases with increasing 
Q 2 . The prediction obtained for Rw is slightly different from that of ref. 
because we are now using a more reliable set of parton distributions. 

4 Quark, antiquark, gluon helicity distribu¬ 
tions and experimental tests 

Since our approach is based on the direct construction of the quark (anti¬ 
quark) parton distributions of a given helicity q± al and we have nothing to 
add to the results of section 2 to obtain A q va i = q^ al — q~ al and A q = q + — q~ 
at Q 2 = Qq. We recall the simple relations used in ref.[|1|] namely 


xAu va i(x,Ql) = xu va i(x,Q 2 0 ) - xd vai (x,Q 2 0 ), (14-a) 

xAd v <a(x,Ql) = ~l/3xd va i(x,Ql), (14-b) 

xAu(x,Ql) = xu(x,Ql) — xd(x,Ql), (14-c) 

xAd(x,Ql) = xAs(x, Ql) = xAs(x, Q%) = 0 . (14-d) 


We show in Fig.13 the different ratios A q/q at Q 2 = 3GeV 2 versus x. One 
finds the usual features for A u/u (positive and growing at large x) and for 
A d/d (negative and growing at large x), whereas A u/u is very large and neg¬ 
ative, which is less conventional, and A d/d is zero. Concerning the helicity 
distribution of the gluon xAG(x,Ql), in the absence of any serious experi¬ 
mental indication, we can only make some speculations. One possibility is 
what one can call the soft gluon , for which one assumes 

xAG(x, Qq) = x 2 G(x, Ql) (15-a) 

and another one, is the hard gluon , where one takes 

xAG(x, Ql) = A(1 — x)x 2 G(x, Ql) (15-b) 


125 1 
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with A = 4 which is the maximum value allowed to obey positivity, i.e. 
AG < G. The integral of AG(x) for the hard gluon case is about four times 
larger than for the soft gluon case. Note that for these two possibilities we 
assume A G(x) > 0, which need not be true. In Fig. 14, together with xG(x), 
we show xAG(x) for these two choices which do not involve any additional 
parameter. 

Having fixed all these hclicity distributions at Q 2 = Qq, at higher Q 2 
they are obtained from the DGLAP evolution equations |5|. For the valence 
contributions, since the splitting functions are the same for q va i and A q va i, 
one can evolve indifferently the r.h.s. or the l.h.s. of eqs. (|14-a|J14-b| ), it will 
lead to the same result which is, anyway, independent of the gluon hclicity 
distribution. This is not the case for the antiquarks and in particular for A u, 
if one evolves the r.h.s. of eq.( |14-c| ), the contribution coming from G(x, Q 2 ) 
will cancel in the difference. This situation corresponds to the assumption 
AG(x, Q 2 ) = 0 for all Q 2 . On the other hand, if one evolves correctly the 
l.h.s. of eq.( |14-c|) , it will depend on AG(x, Q 2 ) for which we can make either 
of the two different choices mentioned above. For the antiquarks, one can 
also assume that eq. (|14-d|) remains valid for all Q 2 or starting from zero at 
Q 2 = Q q one generates a non-zero A q at higher Q 2 from the correct Q 2 
evolution with either the soft or the hard gluon. In Fig. 15a,b we show the 
different ratios A q/q versus x evaluated at Q 2 = M 2 V for the three different 
situations described above. For A u/u the results are very similar and Au/u 
decreases in magnitude for a harder gluon. This is also the case for A d/d, 
but it is reversed for A d/d. Finally we would like to recall that for the gluon 
helicity distribution, due to the QCD evolution, one obtains 
small x region 


2q| in the very 


A G(x,Q 2 )/G(x,Q 2 ) ~ x exp 


0.8 ^JS(Q 2 )lnl/x 


(16) 


for six quark flavors, where S(Q 2 ) = £nt/t 0 and t = (tnQ 2 /A 2 , so at fixed x 
it grows for large Q 2 . 

We now turn to some of the predictions we can make from these dis¬ 
tributions for hclicity dependent observables. Let us start with the polar¬ 
ized structure functions for proton g\[x, Q 2 ), for neutron g™(x, Q 2 ) and for 
deuteron g±(x,Q 2 ). In Fig. 16 we have collected different sets of data for 
xgf, the earlier EMC results [[27j and the SMC results |28j corresponding to 
< Q 2 >= lOGeH 2 and the very recent E143 results |29|] corresponding to 
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< Q 2 >= 3 GeV 2 . We also show for comparison our theoretical predictions, 
down to the lowest x range evaluated at these two Q 2 values. The agreement 
is very satisfactory and gives strong support to our simple relations eqs. (|14-a , 
|14-b| , |14-c| , |14-d|) between unpolarized and polarized parton densities. For the 
integrals of g\ we find 

[ gi(x, Q 2 )dx = 0.132 at Q 2 = lQGeV 2 (17-a) 

J 0.003 


compared to the SMC value [28] 0.131 ± 0.011 ± 0.011 and 


[ gf(x, Q 2 )dx = 0.112 at Q 2 = 3GeV 2 (17-b) 

J 0.029 


compared to the E143 value [29] 0.120 ± 0.004 ± 0.008. In both cases one 


finds results, at least two standard deviations, below the Ellis-Jaffe sum rule 
Idf| with QCD corrections pi] . 


Concerning the neutron polarized structure function xg^(x), we show in 
Fig.17 a comparison of the E142 data |J2] at Q 2 = 2 GeV 2 with our theoret¬ 
ical calculations. The dashed line corresponds to the case where one would 
assume that d quarks are not polarized, i.e. A d(x) = 0, and it clearly dis¬ 
agrees with the data. However by including the d valence quark polarization 
only according to eq. (|14-b|) , we obtain the solid line in perfect agreement 
with the data and we find for Q 2 = 2 GeV 2 


f 1 g™ix)dx = -0.020 . 
Jo 


(18) 


Finally for the deuteron polarized structure function gf, let us recall that 
we have the standard relation 


ai(x) = 2 (sfW + JiW) (1 - 1-5ud) , 


(19) 


where u>d is the H-state probability in the deuteron. We show in Fig.18 
our theoretical prediction at Q 2 = 3 GeV 2 compared to the preliminary E143 
data 11 which are more accurate than the earlier SMC data |34|]. For the 
integral in the x range covered by E143 we find 


(■ 0.8 


' 0.029 


g^(x)dx = 0.043 


( 20 ) 
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in fair agreement with the experimental value [[I3|] 0.044 ± 0.004 ± 0.004. 

Concerning the important issue of the validity of the Bjorken sum rule 
5^], our calculations at Q 2 = 5 GeV 2 leads to 


c i r 


g p i(*,Q 2 




x, Q 5 


dx = 0.158 


( 21 ) 


perfectly compatible with the best experimental estimate [£8| that is 0.166 ± 
0.017. However this is in slight disagreement with the present theoretical 
estimate, including QCD corrections |I6| up to the third order in a s , which 
gives 0.185 ± 0.004. We recall that, to begin with, we did not impose this 
constraint on our distributions. 

Before moving to hadronic collisions, we would like to close this discus¬ 
sion with a few words on the ’’transverse” spin dependent structure function 
92 (x,Q 2 ) which can be measured in polarized deep-inelastic scattering with 
a transversely polarized target. The properties of g 2 have been reviewed re¬ 
cently together with some estimates |57|| and we recall that g 2 can be thought 
as a sum of two terms 


'lit*. Q‘) = Q‘) + 9i(x, Q‘) 


( 22 ) 


where the first term is a twist-2 contribution 
gi(x, Q 2 ) since we have 


|38H determined entirely by 


„ww 

i/2 


(x,Q 2 ) = -gi(x,Q 2 ) + I —gi(y,Q 2 ) 

Jx y 


(23) 


The second term in eq.(22) has twist-3 contributions, determined by quark- 
gluon interactions, and has a priori, no reason to be small. Given our present 
knowledge on g ly we can use eq.(|3D to evaluate g^ w . This is shown in Fig.19 
for the proton case and a comparison with future experimental results will 
provide an estimate of g 2 . 

Next we propose some tests using hadronic collisions, in particular in the 
framework of the future spin physics programme at RHIC which is planned to 
be used as a polarized pp collider |24|] . Many helicity dependent observables 


have been calculated in ref. |^5[ and our purpose here is just to show how to 
update our earlier predictions, using this new set of helicity distributions. 
Let us first consider the parity-violating helicity asymmetry Al defined as 


A i = 


dcr_ /dy — da + / dy 
dcr_/dy + da + /dy 


(24) 
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In W + production, it reads 

A u(x a , Myy)d(x b , Mw) 


MV) = 


(u 


d) 


u(x a , M$y)d(x b , Mw) + {u d) 


(25) 


assuming the proton a is polarized. For W~ production, the quark flavors 
are interchanged. Using our set of polarized quark densities, we find at A /s = 
500GeU, the predictions shown in Fig.20. Near y = +1, A// + ~ A u/u and 
M ~ A d/d evaluated at x — 0.435, whereas near y — — 1 AZ + ~ -A d/d 
and A™ ~ —A u/u evaluated at x = 0.059. Therefore the trends in Fig.20 
can be easily compared with the shapes displayed in Figs. 15a,b which also 
show the sensitivity to the choice of either soft gluon or hard gluon hclicity 
distribution. 

In pp collisions where both proton beams are polarized, there is another 
observable which is very sensitive to antiquark polarizations, that is the 
parity-conserving double hclicity asymmetry A LL defined as 


A ll = 


da ++ /dy + da _/ dy — da + _ / dy — da_ + /dy 

da ++ /dy + da _ /dy + da + _/dy + da/dy 

In lepton-pair production, it reads 


A ll = 


Ei e? [Ag,(x„) Aqi{x b ) + ^ 0,1 

Ei e? fefeOfifo) + « i t )] 


(26) 


(27) 


where e; is the electric charge of the quark qy 

We have calculated All at a/s = lOOGeU which seems best appropriate 
to the acceptance of the detectors at RHIC and the results, in the soft gluon 
case, are shown in Fig.21. We observe that All is positive because it is 
dominated by Aw Ah, where A u is negative. All increases with increasing 
lepton-pair M and of course for A u = Ad = 0 we would have A L l = 0. 
Obviously these tests can be extended by many other examples as shown in 
ref. [|25| , in particular parity-violating and parity-conserving hclicity asymme¬ 
tries in W ± , Z production in pp and pn collisions. 


5 Double spin transverse asymmetries Aj'j' 

So far we have considered collisions involving only longitudinally polarized 
proton beams, but of course at RHIC, transversely polarized protons will 
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be available as well Bl. This new possibility is extremely appealing be¬ 


cause of recent progress in understanding transverse spin effects in QCD, 
both at leading twist |3D|] and higher twist levels m • For the case of the 
nucleon’s helicity, its distribution among the various quarks and antiquarks 
can be obtained in polarized deep-inelastic scattering from the measurement 
of the structure function g\(x) mentioned above. However this is not pos¬ 
sible for the transversity distribution h\(x) which describes the state of a 
quark (antiquark) in a transversely polarized nucleon. The reason is that 
hi(x), which measures the correlation between right-handed and left-handed 
quarks, decouples from deep-inelastic scattering. Indeed like gi(x), h\(x) is 
leading - twist and it can be measured in Drell-Yan lepton-pair production 
with both initial proton beams transversely polarized pi. Other possibilities 


have been suggested ^ 1 | but in the framework of this paper, we will envisage 
also a practical way to determine hi(x), by using gauge boson production in 
pp collisions with protons transversely polarized. Let us consider the double 
spin transverse asymmetry defined as 


Att — 


^tt ~ a n 

°tt + a n 


(28) 


where crtT(^Ti) denotes the cross section with the two initial protons trans¬ 
versely polarized in the same (opposite) direction. Assuming that the under¬ 
lying parton subprocess is quark-antiquark annihilation, we easily find for Z 
production 


Att — 


Ei=u,d (X 2 - at) \hT(x a )hf(x b ) + (x a X b ) 


J2i=u,d K 2 + ti-) [ qi(x a )qi(x b ) + (x a <-»• x b )] 


(29) 


This result generalizes the case of the lepton-pair production |39|] through 
an off-shell photon 7 * corresponding to b{ = 0 and cq = e* and which gives 


Att — cltt~ 


Ei \hf(x a )hf (x b ) + (x a ^ x b ) 


Ei e f [qi(x a )qi(x b ) + (x a <-> x b )} 


(30) 


where Ott is the parton asymmetry which has a simple expression in the c.m. 
frame of the lepton-pair j[k]| . 

For W ± production, which is pure left-handed and therefore does not 
allow right-left interference, we expect Att = 0 , since in this case af = bj. 
This result is worth checking experimentally. 
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So far there is no experimental data on these distributions h\(x) (or 
h\(x)), but there are some attempts to calculate them either in the frame¬ 
work of the MIT bag model [[IIJ or by means of QCD sum rules [^]. However 
the use of positivity yields to derive a model-independent constraint on h\ (x) 
which restricts substantially the domain of allowed values |43|]. Indeed one 
has obtained 

q(x) + A q[x) > 2\h\(x)\ . (31) 

which is much less trivial than 


q(x) > |fc?(l)l 


(32) 


as proposed earlier in ref. [39 . 


In the MIT bag model, let us recall that these distributions read [39 


q = f + g 2 , 


Aq = f - 1/3 g 2 and h\ = f + l/3 ^ 2 


(33) 


and they saturate eq.(pl]). In this case, we observe that h\{x) > A q(x) but 
this situation cannot be very general because of eq. (0)- As an example let us 
assume h\(x) = 2A q(x). Such a relation cannot hold for all x and we see that 
eq.(|3Tp, in particular if A q(x) > 0, implies q(x) > 3A q(x). This is certainly 
not satisfied for all x by the present determination of the u quark hclicity 
distribution, in particular for large x where Af(x) is large p7|, [28], |29|] . The 
simplifying assumption h\(x) = A q(x), based on the non-relativistic quark 
model, which has been used in some recent calculations PH PT| is also not 


acceptable for all x values if A q(x) < 0 because of eq. (|3T| ). To illustrate 
the practical use of eq.(^l]), let us consider eqs. (|14-a|) and (|14-d) . It is then 
possible to obtain the allowed range of values for h“(x), namely 


u(x) - -d(x) > I h1(x)\ 


(34) 


and similarly for u. 

We now turn to some predictions for A TT which will depend heavily on 
our assumptions. Since the u quark term is expected to dominate because of 
the charge factor, we assume /if = 0 and we take the equality sign in eq.(34) 
assuming /if(x) > 0 and /rf(x) < 0. We show in Figs.22a,b our results for 
Att/o,tt in dilepton production at two energies y/s = 50 and 100 GeV for 
several values of the dilepton mass M. The effect is larger for increasing M 


15 









































because at fixed energy it corresponds to higher x values where A u and A u 
are larger, as well as h™ and h Of course, at fixed M the effect decreases 
with increasing energy. Finally we show in Fig.23 our predictions for A TT 
in Z production at two different energies = 350 and 500 GeV. Clearly 
these predictions are only a guide for future experiments at RHIC which will 
indeed lead to the actual determination of h\{x). 


6 Concluding remarks 

We have presented a new set of quark-parton densities in terms of Ferrni- 
Dirac distributions depending on nine free parameters which were determined 
from the spin-average structure functions at low Q 2 . Some simple relations 
between unpolarized and polarized densities, which were postulated earlier, 
are well supported by recent data on spin-dependent structure functions. We 
have also proposed a simple expression for the gluon density in terms of a 
Bose-Einstein distribution with no additional free parameter. We have then 
used a straightforward DGLAP Q 2 evolution to get access to a broad x and 
Q 2 range in order to test our approach with a large number of deep-inelastic 
scattering data. The predictions give a very satisfactory description of the 
CCFR, NMC, BCDMS, SLAC data and we also get, in a natural way, the 
sharp rise of F.^^x, Q 2 ) for small x, recently observed at HERA by the HI 
and Zeus Collaborations. So this is a brillant confirmation of the validity of 
perturbative QCD, to be confronted to more accurate future data at HERA. 
We have also discussed the relevant question of flavor symmetry breaking of 
the light sea quarks (or antiquarks) of the nucleon which has to be measured 
in a broader x and Q 2 range. To extend the polarized antiquark (or sea quark) 
parton densities to large Q 2 , we have to make a choice for the polarized gluon 
distribution and we used two possiblities, the soft gluon and the hard gluon. 
None of these is favoured at present and experiment will have to decide. 
Predictions were given for future measurements in hadronic collisions with 
longitudinally and transversely polarized proton beams at RHIC which will 
allow further tests of our approach in the spin sector, in particular for dilepton 
and W ± , Z production. 
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Figure Captions 

Fig.la The structure function xF^ N (x, Q 2 ) versus x. Data are from ref.0 
at Q 2 = 3 GeV 2 , the solid line is the result of our fit and the dashed 
curve represents the MRS (A) parametrization at Q 2 = 4GeV 2 from 
ref. 0 . 

Fig.lb The antiquark contribution xq(x ) = xu(x) + xd(x) + xs(x) at Q 2 = 
3 GeV 2 (full circles) and Q 2 = 5GeV 2 (full triangles) versus x. Data 
are from ref. 0 , the solid line is the result of our fit and the dashed 
curve represents the MRS (A) parametrization at Q 2 = 4 GeV 2 from 
ref. 0 . 

Fig.2a The difference Tf (x)—F^(x) at Q 2 = 4GeU 2 versus x. Data are from 
ref. |§], the solid line is the result of our calculation and the dashed 
curve represents the MRS (A) parametrization at Q 2 = 4 GeV 2 from 
ref. H . 

Fig.2b The ratio F2(x)/F$(x) at Q 2 = 4 GeV 2 versus x. Data are from 
ref.||, the solid line is the result of our calculation and the dashed 
curve represents the MRS (A) parametrization at Q 2 = 4GeV 2 from 
ref. U . 

Fig.3 The unpolarized parton distributions as a function of x at two differ¬ 
ent Q 2 values, Q 2 = 3 GeV 2 (solid lines) and Q 2 = 20 GeV 2 (dashed 
lines). 

Fig.4 The gluon density xG(x, Q 2 ) versus x at two different Q 2 values, 
Q 2 = 3 GeV 2 (solid line) and Q 2 = 20GeU 2 (dashed line). 
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Fig.5 


Fig.6a 


Fig.6b 
Fig.7a 


Fig.7b 
Fig.7c 
Fig.8a 


Fig.8b 
Fig.9 

Fig.10 

Fig.11 

Fig.12 
Fig.13 


The description of the CCFR data [Q of the xF^ N (x,Q 2 ) structure 
function, given by our parton densities. The theoretical curves are 
shown after an overall renormalization by a factor 1.04. 


The description of the NMC data |T(J (full circles) and BCDMS 
data (U| (open circles) of the F/ P (a;, Q 2 ) structure function given 
by our parton densities versus Q 2 for x bins between x = 0.0125 
and x = 0.14. The theoretical curves are shown after an overall 
renormalization by a factor 1.12. 


Same as Fig.6a for x bins between x = 0.18 and x = 0.75. 


The description of the Zeus data |T^| (full triangles) and HI data |T^ 
(full circles) of the F- 2 P (x, Q 2 ) structure function given by our parton 
densities versus x for Q 2 bins between Q 2 = 4.5 GeV 2 and Q 2 = 
25 GeV 2 . 


Same as Fig. 7a for Q 2 bins between Q 2 = 35 GeV 2 and 200GeR 2 . 
Same as Fig.7a for Q 2 bins between Q 2 = 250GeR 2 and 5000GeR 2 . 


The description of the NMC data fllj (full circles) and SLAC data 
[0 ''open circles) of the Q 2 ) structure function in the low 

Q 2 , low x region given by our parton densities for x bins between 
x = 0.00375 and x = 0.0175. The theoretical curves are shown after 
an overall renormalization by a factor 1.12. 


Same as Fig.8a for x bins between x = 0.025 and x = 0.14. 


The W + rapidity distribution in pp collisions at a/s = 1.8TeR pre¬ 
dicted by our parton densities. 


Our theoretical prediction for A^y 
compared to the NA51 data |22| . 


(see eq.(|T0|)) versus y/r = M/a/s 


Our predicted ratio (d — u)/(d + u ) versus x for two different Q 2 
values ; solid line Q 2 = 25 GeV 2 and dashed line Q 2 = M 


The ratio Rw (see eq.(p^)) versus y at y/s = 250 and 350 GeV. 

The ratios A u/u, A d/d, A u/u and A d/d versus x at Q 2 = 3 GeV 2 . 
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Fig.14 

Fig.15a 

Fig.15b 
Fig.16 

Fig.17 

Fig.18 
Fig.19 
Fig.20 

Fig.21 

Fig.22a 

Fig.22b 
Fig.23 


The hclicity gluon distribution xAG(x) (soft gluon, dashed line and 
hard gluon, dotted line) versus x at Q 2 = 3 GeV 2 and for comparison 
the unpolarized gluon distribution xG(x) (solid line). 

The predicted ratios A u/u and A u/u versus x at Q 2 = M(y. The 
three curves correspond to the three cases for AG: identically zero 
(solid line), soft gluon (eq.( |15-a|) ) (dashed-dotted line) and hard 
gluon (eq. (|15-b|) ) (dashed line). 

Same as Fig. 15a for A d/d and A d/d. 

xg\(x ) versus x at < Q 2 >= 3 GeV 2 and 10 GeV 2 . Data are from 
ref.|27j (full squares), ref. |28|] (full circles) and ref. |29j (open circles). 
The curves are theoretical predictions at Q 2 = 3 GeV 2 (dashed line) 
and Q 2 = 10 GeV 2 (solid line). 


xg™ (x) versus x at < Q 2 >= 2 GeV 2 . Data are from ref. |32j together 
with our theoretical predictions (dashed line is the contribution of 
A u(x) and A u{x) only and solid line is the contribution of A«(x), 
A u(x) and A d va i(x)). 

xgf(x) versus x at < Q 2 >= 3 GeV 2 . Prelimining data are from 
ref. El together with our theoretical prediction. 


Our theoretical prediction for g 2 (x,Q 2 ) at Q 2 = 3 GeV 2 for proton 


target assuming g 2 = 0 (see eq. (|22|) ). 

Parity-violating hclicity asymmetry Al versus y for W ± production 
in pp collisions at yfs = 500GeD. Solid lines correspond to soft 
gluon and dashed lines to hard gluon. 


Parity-conserving hclicity asymmetry A^i versus y for dilepton pro¬ 
duction at y/s = lOOGeD and different values of the lepton-pair 
mass. The curves correspond to the soft gluon case. 

Predictions for A tt /citt in dilepton production versus y for y/s = 
50 GeV and M = 5, 10, 15 GeV. The curves correspond to the soft 
gluon case. 

Same as Fig.22a for y/s = 100GeD. 

Predictions for Att in Z production versus y for yfs = 350 and 
500GeD. The curves correspond to the soft gluon case. 
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